Spontaneous motion of a droplet coupled with a chemical wave 
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We propose a novel framework for the spontaneous motion of a droplet coupled with internal 
dynamic patterns generated in a reaction-diffusion system. The spatio-temporal order of the chem- 
ical reaction gives rise to inhomogeneous surface tension and results in self-propulsion driven by the 
surrounding flow due to the Marangoni effect. Numerical calculations of internal patterns together 
with theoretical results of the flow fields at low Reynolds number well reproduces the experimental 
results obtained using a droplet of Belousov-Zhabotinsky (BZ) reaction medium. 
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Spatio-temporal patterns are widely seen in living sys- 
tems; target, spiral, stripe, and dot patterns have been 
observed at various scales from the interior of a cell to 
a swarm of cells. Most of studies have been focused on 
patterns at larger scales, which can be successfully re- 
produced using reaction-diffusion dynamics [lj. In con- 
trast, it is only recently that internal patterns in a single 
cell have been visualized. These patterns are expected 
to relate to cellular functions; examples include calcium 
ions for signal transduction [2j], Min proteins for cell di- 
vision Q, and actin cytoskelton for mechanical proper- 
ties [H . Although pattern formation in a cell is expected 
to be analyzed in the framework of a reaction-diffusion 
system, as demonstrated in in vitro experiments |3(, suf- 
ficient understanding on the connection between pattern 
formation and cellular function is awaited. In this paper, 
we focus motility, as a typical aspect of cellular functions, 
arising from internal patterns. 

Several artificial systems imitating cell motility have 
been proposed as self-propelled particles |5(. Although 
no external force is exerted on the particles (force-free 
condition) , the motion is induced by the asymmetric dis- 
tribution of an electric field, concentration of chemicals, 
temperature, and so on. These asymmetric distributions 
are either a priori embedded in the asymmetry of the sur- 
face properties of the self-propelled objects [6J , or a pos- 
teriori created by nonlinear effects; motion itself destabi- 
lizes symmetric distribution, for instance, through advec- 
tive flow 0. In both cases, however, most studies have 
focused on motion under steady distributions. 

In order to understand the dynamic features of cell 
motility, a system connecting dynamic pattern with mo- 
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tion is desirable. In fact, experimental and numerical 
evidence of chemo-mechanical coupling in such systems 
has been demonstrated [8|, [9j . In this letter, we propose 
a theoretical framework for a chemical system exhibit- 
ing self-organized patterns, leading to spontaneous mo- 
tion. We consider that the Marangoni effect is suitable 
for this purpose as it has been shown to drive an object 
under force-free conditions by an inhomogeneous interfa- 
cial tension arising from a gradient in chemical concen- 
tration 11 1 . In our system, the energy supply and 
consumption can generate a pattern in a droplet through 
nonlinear chemical kinetics, and the pattern at the in- 
terface of the droplet creates inhomogeneous interfacial 
tension. This generates a flow surrounding the droplet, 
resulting in motion. 

We discuss a spherical droplet of incompressible fluid 
with a radius of R in another fluid. We consider the 
Stokes equation under the approximation of low Reynolds 
number [12j. The inertia term pdv/dt is also neglected 
where p is fluid density. We will justify these assump- 
tions later using experimental values. The concentra- 
tion of chemical species, which determines the interfacial 
tension, obeys the reaction-diffusion-advection equation. 
The sets of governing equations are described as 

pd t v = - V> + ?/V 2 v = 0, (1) 
V ■ v = 0, (2) 
d t c + (v • V) c = F (c) + V (L»Vc) , (3) 

where v, p, and r\ are the flow velocity, pressure, and vis- 
cosity of the fluid, respectively. c(r) is a vector composed 
of concentrations of chemicals, F(c) corresponds to the 
reaction kinetics, and D is a diffusion coefficient. 

The flow velocity is coupled with the chemical reaction 
through interfacial tension at the droplet interface. This 
coupling is described as a force balance in a tangential 



direction at an interface [111 Il3j : a. 
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(dj/d6)/R, where a is the stress tensor, and j(9) is the 
interfacial tension profile, which is dependent on c. 
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FIG. 1: (a) Schematic representation of the system under 
consideration. The broken line corresponds to the symmetry 
axis, (b) Setting of the coordinates in the droplet system. 



For simplicity, we only consider the axisymmetric sys- 
tem schematically shown in Fig. HJa). Since the bound- 
ary condition is given at the droplet interface, it is conve- 
nient to set the spherical coordinates so that the droplet 
is fixed, as shown in Fig. [T^b). We have to adopt the 
boundary condition that Vr — w = and Vg = vi°^ 
at r = R and v — > vo = u as r oo, where r is the 
distance from the droplet center, and u is the droplet 
velocity in a laboratory system. The superscripts "(i)" 
and "(o)" correspond to the fluids inside and outside of 
the droplet, respectively. From the symmetric property, 
we aligned u with positive z; i.e., u = ue z . We also have 
to consider the condition that the force exerted on the 
droplet, f, is zero; i.e., f = /d5n • a-^°'\ r _ R = 0. 

By solving Eqs. ((T|) and © with the above conditions, 
the droplet velocity is calculated as 



9t?W + 6t7(°) 



-aT 1 , 



(4) 



where Ti is the first-mode coefficient of the Legendre 
expansion of 7 (0) 0Q: 7 (0) = £~ =0 r„P„(cos 9), 
where P n is the Legendre poly nominal of order n. The 
higher modes (n > 2) contribute only to deformation 
and are independent to the translational motion for small 
deformation. The detailed derivation on the solutions of 
the Stokes equation is shown in Supplemental Materials 
(SM) [3. 

We focus on a specific system using the Belousov- 
Zhabotinsky (BZ) reaction, which is common in exper- 
imental systems that can exhibit spatio-temporal pat- 
tern formation, such as target and spiral patterns in a 
two-dimensional system (l6j and a scroll ring in a three- 
dimensional system 17] . As a mathematical model for 
the BZ reaction, we adopt the Oregonator [l8j], which 
is widely used due to its reliability as demonstrated in 
physico-chemical discussion on the elemental processes 
of chemical reactions. In the Oregonator, c is composed 
of the two variables U and V, which corrrepond to the 
concentrations of HBrC>2 and oxidized catalyst, also re- 
ferred to as the activator and inhibitor, respectively. F(c) 




FIG. 2: (Color online) Numerical calculations for a single 
pulse, (a) Snapshots of V at the cross section and displace- 
ment of a droplet. Blue (light gray) and red (dark gray) colors 
correspond to higher and lower V, respectively. The displace- 
ment is shown as the vertical position. At t = 0, a chemical 
wave was initiated by setting U — 1 at a certain region in the 
lower half of the droplet, (b) r n (1 < n < 4) against time, 
(c) Velocity u and the position z of the droplet against time. 
The parameters: R = 125, D = 1600. 



is described as: 



[7(1 -U)-fV 

u-v 



U-q 
U + q 



(5) 



where e, q, and / are the parameters that determine the 
characterisitcs of the BZ reaction. The Oregonator model 
is nondimensionalized, so that the time unit in the cal- 
culation is set as T. 

The reaction-diffusion- advection equation (Eq. ([3])) is 
numerically solved inside a droplet, which corresponds to 
the experimental conditions in which a droplet of BZ re- 
action medium (BZ droplet) is inside another fluid. Here, 
by taking the axisymmetry into consideration, the calcu- 
lations can be performed in a two-dimensional field . We 
adopt a monotonous increasing function for the interfa- 
cial tension against V: 



j(V) = 7o + kV. 



(6) 



This assumption is supported by the experimental mea- 
surements; the air/water interfacial tension in the ferroin- 
catalyzed BZ reaction medium is higher in the oxidized 
state than that in the reduced state [l^, [2(| . In addition, 
Marangoni flow was observed at the oil/ water interface, 
as well as the air/water interface j^j. Note that the coef- 
ficient may include the effect due to the heat generation 
by reaction, which only reduces its absolute value. 

The parameters for the Oregonator model were set as 
q = 0.001, e = 0.05, and / = 2.5, which correspond to 
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FIG. 3: (Color online) Snapshots of streamlines and stream 
function at the cross section corresponding to the snapshots in 
Fig. [2] (a). Red (dark gray) and blue (light gray) colors corre- 
spond to positive and negative values of the stream function, 
respectively [lq] . The direction of the rolls in each droplet is 
shown by red (dark gray) or blue (light gray) coloring (clock- 
wise: red (dark gray) in the left hemisphere and blue (light 
gray) in the right one, counterclockwise: blue (light gray) in 
the left hemisphere and red (dark gray) in the right one). 



the excitable condition, in which a chemical wave prop- 
agates only from the initiated point. In the numerical 
calculation, the velocity field was calculated by summing 
the modes of n = 1 to 8. We confirmed that the cut-off 
modes did not critically affect the numerical results. 

We performed numerical calculations under two typical 
sets of conditions; one is a single pulse that corresponds 
to a target pattern, and the other is a scroll ring 



17]. 



which is a natural extension of a spiral wave in a three- 
dimensional medium. 

Figure [2] shows the results of numerical calculations for 
a single pulse. At t = 0, a chemical wave was initiated in 
the lower half of the droplet, and then a single chemical 
wave propagated isotropically. Figure [5] (a) shows snap- 
shots of V in a cross section of the droplet, and Fig. [2] (b) 
and (c) are the time series of T n (1 < n < 4), velocity and 
position of the droplet, respectively. The chemical wave 
arose from the lower part of the droplet, and propagated 
in a circular shape. When the wave reached the lower 
boundary of the droplet, the droplet began to move up- 
ward. When the chemical wave reached the upper bound- 
ary, the droplet moved back in the opposite direction. 
Nevertheless, there was a net motion in a cycle (see Fig. 
[DJc) and movies in SM (lBj.) Figure |3] shows snapshots 
of streamlines in a cross section inside a droplet from the 
droplet-fixed frame. The direction of the convective rolls 
was inverted at the moment the droplet moved back. 

In the numerical calculation for a scroll ring, the pat- 
tern was prepared using the following standard proce- 
dure; a chemical wave was first initiated, and a part of 
the chemical wave was omitted by setting U and V as 
steady-state values at t ~ 2.5T. The numerical results 
for the snapshots and time series of the position of the 
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FIG. 4: (Color online) Numerical calculations for a scroll 
wave, (a) Snapshots of V at the cross section and displace- 
ment of a droplet, shown in the same manner as in Fig. [2] 
At t = 0, a chemical wave was initiated by setting U = 1 
at a certain region in the lower half of the droplet, and the 
chemical wave at the central region of the droplet was omitted 
at t — 2.5T. (b) Time series of the position of the droplet. 
Parameters: R = 125, D = 60. 



droplet are shown in Fig. 01 The droplet moved back 
and forth, gradually moving upward, coupling with the 
development of the scroll ring (see the movie in SM [15|). 
In Fig. H|a), the chemical wave first reached the lower 
half of the droplet, which lead to upward motion. As 
the point of contact moved toward the upper half of the 
droplet, downward motion resulted. The droplet finally 
stopped, since the filament (alignment of the phase singu- 
larity points) shrank and disappeared due to its intrinsic 
instability [ill ]. 

In order to compare the theoretical and experimental 
results, we undertook experiments on spontaneous mo- 
tion of a BZ droplet inside an oil phase as schematically 
shown in Fig. 03a). We put a droplet of BZ medium 
with a volume of 1 /il into an oil phase (oleic acid), in 
a petri dish made of polytetrafluoroethylene to prevent 
the droplet from coming into contact with the bottom of 
the petri dish. We then observed the BZ droplet from 
above using a digital video camera (DCR-HC62, Sony, 
Japan) equipped with a close-up lens (CM-3500, Raynox, 
Japan). The composition of the BZ reaction medium was 
[NaBr0 3 ] = 0.3 M, [H 2 S0 4 ] = 0.6M, [CH 2 (COOH) 2 ] = 
0.1M, [NaBr] = 30 mM and [Fe(phen) 3 S0 4 ] = 5mM. 

The experimental results are summarized in Fig. [SJb) 
and (c). The BZ droplet moves in a horizontal plane. 
The axis of the position in (b) is set so that the upper 
direction in (c) is positive. As the BZ reaction medium 
was prepared in an oscillatory condition, a chemical wave 
was initiated spontaneously from the point determined 
stochastically. In Fig. [5jc), a chemical wave was initiated 
in the lower half of the droplet. When the propagating 
chemical wave touched the lower interface, the droplet 
began to move upwards. It then moved back in the op- 
posite direction when the chemical wave propagated over 
the entire droplet. 
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FIG. 5: (Color online) Experiments on the spontaneous mo- 
tion of a BZ droplet, (a) Schematic representation of the 
experimental setup, (b) Time change in the center of mass 
of the BZ droplet, (c) Snapshots in a horizontal plane (x-z 
plane) at every 1 s. Red (dark gray) and blue (light gray) 
regions correspond to the oxidized and reduced states in the 
BZ reaction, respectively. 



From the experimental results, the maximum droplet 
velocity was measured as u max = 0.1 mms -1 . The viscos- 
ity of the BZ medium, rjW, is ~ 10 -3 kgm^ 1 s^ 1 and that 



The Reynolds 



of oleic acid, T]^°>, is ~ 10~ 2 kg uT 
number in the present experimental system is calculated 
with the kinetic viscosity, v, as Re = uR/v ~ 0.1, 
which shows that the Stokes approximation is applicable. 
The relaxation of the velocity field is characterized by 
t u = R 2 jv ~ Is, while the time scale of the wave prop- 
agation is r p = Rj ' y/D/e ~ Is. With these values, we 
confirm that the time-derivative in Eq. ([T]) does not qual- 
itatively affect our results, but modifies Eq. (j4} within 
15% following the method in 22j- Therefore, to the sake 
of simplicity, we assume the Sow both inside and outside 
of the droplet rapidly changes following the changes in 
chemical concentration. Using these values and the re- 
sults of the numerical calculations, the maximum differ- 
ence in the interfacial tension between the oxidized and 
reduced conditions is estimated as A7 ~ O.OlmNm -1 . 
Although it is difficult to directly measure the interfacial 
tension at the oil/ water interface in this particular sys- 
tem, A7 at an air/water interface of the BZ medium is 
available The value in the reference was ten times 
higher than our estimation for A7 at an oil/ water inter- 
face. In general, the interfacial tension of an oil/water 
interface is lower than that of an air /water interface [23[ , 
and thus the present estimation seems reasonable. 

In summary, we studied droplet motion coupled with 
internal dynamic pattern formation through an interfa- 
cial tension gradient. The Stokes equation was adopted 
for a spherical droplet, and the motion of the droplet 
was calculated, which well reproduced the experimental 
results. Although we have not so far experimentally suc- 
ceeded in generating scrolling rings in a BZ droplet, we 
expect this motion will be realized by using photosensi- 
tive BZ reaction. A similar model could also be applied 
to other types of motion; for example, rotational and 



oscillatory motion. It could further be extended to the 
motion of a nonspherical droplet. 
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